Antilambda production in Au+Au collisions at 11.7 A-GeV/c 
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We present results from Experiment E917 for antilambda and antiproton production in Au+Au 
collisions at 11.7 A-GeV /c. We have measured invariant spectra and yields for both species in central 



and peripheral collisions. We find that the A/p ratio near mid-rapidity increases from 0.26 



+0.19 



peripheral collisions to 3.6_ 
theoretical estimates. 



in central collisions, a value that is substantially larger than current 



PACS numbers: 25.75.-q, 13.85.Ni, 21.65. +f 
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Enhanced antimatter and strangeness production have 
both been proposed as signatures of the phase transition 
from normal hadronic matter to the quark-gluon plasma 
(QGP) |l], ^, E|. The yields of strange antibaryons com- 
bine both of these signatures in a single production chan- 
nel, and may prove to be a particularly sensitive probe 
for the formation of the QGP. Of particular interest is 
the ratio A/p 0, which reflects the relative abundance 
of s quarks to light antiquarks. 

However, both A and p are also produced when the 
collision does not form a QGP. Hadronic matter with 
sufficient rescattering can reach a chemical equilibrium 
containing the full spectrum of hadrons. Equilibrium 
chemical models of heavy- ion collisions are generally able 
to reproduce the measured ratios of total yields of par- 
ticles such as pions and kaons. Such models also predict 
values for the ratio A/p in the range 0.15-0.9 at AGS 
beam energies j|, ||. A thermal model analysis that ac- 
counts for the measured K+/-7T+ and K~/K + ratios || 
gives a most probable value of A/p — 1, with unrealistic 
values yielding an upper limit of 2 for Au+Au collisions. 

Hadronic cascade models which explicitly follow the 
reactions and trajectories of hadrons in the collision zone 
^| predict values for the A/p ratio of less than 1. 
The absorption cross-section for A-N processes is not well 
known at these energies, opening the possibility for dif- 
ferential absorption of the two species to play a role in 
determining the value of the ratio A/p and its rapidity 
dependence. However, a study has shown that a reason- 
able interpretation of this value cannot bring the predic- 



tions of cascade models above a ratio of ~ 1 H . Recent 
data from E877 on the anti-flow of p M are successfully 
described by the hadronic cascades that use the normal 
antiproton annihilation cross-section. 

Enhanced antibaryon production may also be due to 
many-body collisions between pions, e.g. mr — > p + p 
JTo| , pd| . The rate of antibaryon production can be es- 
timated from the inverse reaction p + p — > nir. Cal- 
culations by Rapp and Greiner Jll[ indicate that a 
significant number of both p and A may come from these 
processes, but it is not yet clear how to incorporate these 
many-body collisions into current hadronic cascade cal- 
culations. 

In summary, hadronic models with reasonable inputs 
generally predict A/p < 1. Thus the experimental detec- 
tion of a A/p ratio significantly larger than 1 may signal 
the breakdown of these hadronic models, and possibly the 
onset of other degrees of freedom in the system. Of par- 
ticular interest is the centrality dependence of the A/p 
ratio, which may discriminate Jl2| between the various 
models attempting to describe antibaryon production: 
QGP, hadronic scattering, and thermal models. Canon- 
ical thermal models should be used in this comparison, 
since the volume of the emitting system has been shown 
to affect the ratio of strange to non-strange particles by 
up to a factor of four (b| . 

There are indications of an abnormally large value of 
the A/p ratio from other experiments. Experiment E859 
reported a directly measured A/p ratio of 2.9±0.9(stat)± 
0.5(sys) near mid-rapidity in central Si+Au collisions, at 
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a beam energy of 14.6 A-GeV/c E3] . An indirect estimate 
of A/p > 2.3 (98% C.L.) at y=1.6 for p ± ~ in central 
Au+Pb collisions was obtained at a beam energy of 11.5 
A-GeV / c by comparing inclusive antiproton yields from 
two different experients (E864 and E878) with largely 
different acceptances |l5fl . 

In this letter, we present a direct measurement of A in 
Au+Au collisions near mid-rapidity and over an extended 
range in p± . We also measure the total p spectrum, which 
includes contributions from both direct p production and 
A decay. Using the A measurement, we can estimate the 
contribution to p from A decay and thus extract the yield 
of directly produced antiprotons (pdirect)- We find the 
ratio A/pdirect to be large in central collisions, and much 
lower in more peripheral collisions. 

The experiment was carried out using the magnetic 
tracking spectrometer previously employed in AGS ex- 
periments E802, E859, and E866 as described in detail 
in Refs. |l(], [l7], [l^, |l9). An efficient second- level trigger 
was used in order to study the production of rare parti- 
cles such as antiprotons. A 197 Au beam with momentum 
11.7 A-GeV/c was incident on a 1-g/cm 2 thick Au target. 
The movable spectrometer subtended polar angles from 
19° to 34°. Magnetic field settings of ±0.4 T were used, 
for which the momentum resolution of the spectrometer, 
Sp/p, ranges from 1% to 2%. At these angle and field 
settings, we have acceptance for p and A in the rapid- 
ity range 1.0 < y < 1.4, where mid-rapidity is y = 1.6 
at this beam energy. Particle identification is based on 
the measured momentum combined with a time-of-flight 
(TOF) wall with a resolution of 135 ps situated 6 m from 
the target. 

We define centrality classes based on the fraction of 
the total interaction cross-section (a to t = 6.8 barns pO[). 
Using the total kinetic energy of spectator particles mea- 
sured in a zero-degree calorimeter, the data were divided 
into two centrality classes, called "Central" and "Periph- 
eral", containing 0% - 12% and 12% - 77% of atot, re- 
spectively. We also report a "Minimum Bias" data set 
without centrality cuts. 

The reconstruction, particle identification and tracking 
algorithms have been described in detail elsewhere [ fl8[ . 
Antiprotons are accepted up to momentum 2.9 GeV/c, 
above which the 3a identification bands of the K~ and p 
time of flight spectrum begin to overlap. A lower momen- 
tum cutoff of 0.5 GeV/c was necessitated by hadronic in- 
teraction rates in the spectrometer material. Pions used 
to reconstruct A's were required to have momenta be- 
tween 0.35-1.8 GeV/c. 

A firm understanding of the background in the TOF- 
identified p sample is required for the measurement of 
the p spectrum. Spectra of the time-of-flight residuals 
At = (t rneas - t expected ) I a were created in order to es- 
timate these backgrounds in the p sample. Here t meas 
and t eX p ec t e d are the measured and expected flight times 
of p's to the TOF wall, and a is a particle-by-particle 



estimate of the timing resolution based on the TOF wall 
resolution combined with a Monte Carlo simulation of the 
tracking in the spectrometer. The At distributions are 
parameterized by a Gaussian for the true p signal and 
a background with two components: an exponentially 
falling contribution from the tails of K - and 7r~ bands, 
and a flat random contribution due to misidentified parti- 
cles. The backgrounds extracted range from 15% to 75% 
of the total signal, with the largest backgrounds at the 
extremes of the momentum ranges. Rapidity-transverse 
mass bins for which the background fraction exceeds 50% 
have been excluded from further analysis. 

The counting statistics in determining the background 
are included in the reported statistical error bars. We 
estimate an additional systematic error in the extracted 
p yield due to background corrections and acceptance 
uncertainties to be 10%, and the error in the p inverse 
slope, arising mostly from acceptance uncertainties, to 
be 5%. The measured p invariant spectra are shown as 
open circles in Fig. | as a function of transverse mass 
(mj_ = \/p\ + tuq) , along with independent fits (dashed 
lines) of a Boltzmann form, B(T), i.e. 

1 d 2 N = dN 
2iTm± dm±dy dy 

dN 1 exp(— (mj_ — too)/T) 
~dy~2^ m± m\T + 2m T 2 + 2T 3 ' 

The parameters of the fit to the total p spectra 
(dNp total j dy , Tp total ) are given in columns two and three 
of Table |. These spectra include contributions from 
both direct production (pdirect) and anti-hyperon decay 
{p decay), which we now consider. 

We reconstruct A's through direct measurement of 
(p, 7r + ) pairs. This decay channel has a branching ratio 
of 63.9% |plf . The spectrum of reconstructed invariant 
mass for accepted (p, ir + ) pairs is shown in Fig. |2[ The 
experimental resolution for reconstructing A's is found 
by initially reconstructing the A invariant mass spec- 
trum in the same data set. The peak of the A spec- 
trum is found to beat 1.1161±0.0004 GeV/c 2 , consistent 
with the accepted value of ttia = 1.1157 GeV/c 2 . The 
Gaussian width of the peak is a — 1.3 MeV/c 2 , which 
is attributable to experimental resolution. To define A 
candidates, a cut of ±3<r is taken around the nominal 
m\. To estimate the background under the A peak, we 
have constructed a mixed-event background, shown as 
the solid histogram in Fig. |^. Pairs are constructed by 
selecting random p's and ir + 's from the set of all pairs 
in the centrality class being considered. The background 
is normalized to the region outside the A peak, defined 
as more than 6<r from m\. Both real and mixed pairs 
were included only if their opening angle was greater than 
15 mrad. This removes most of the inefficiency for de- 
tecting close tracks. The remaining two-track inefficiency 
is estimated by comparison with the mixed-event sample. 
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The correction is then applied on a pair-by-pair basis in 
generating the background invariant mass spectrum as 
well as to the signal's final m± spectra. This correction 
has less than a 5% effect on the yields. 

We further correct the estimated background for the 
residual correlations arising from the presence of particles 
from signal pairs in the mixed-event background. Pairs 
in the mixed-event background with a member that came 
from the signal region are assigned a weight equal to the 
background fraction. This must be done in an iterative 
fashion; four iterations are sufficient to attain conver- 
gence. This correction increases the extracted A yields 
by approximately 15%. From the small remaining dis- 
crepancy between the shape of the mixed-event sample 
and the actual sample, we estimate the systematic error 
on the measured A yield arising from background sub- 
traction to be 15%. The A invariant spectra are shown 
as solid circles in Fig. [I] as a function of transverse mass 
in all centrality classes. 

Since the p spectra have a strong contribution from 
A decays, it is advantageous to perform a simultaneous 
fit to both spectra, in order to obtain the most accu- 
rate measure of the A/p ratio. We have used a Monte 
Carlo study to parameterize the distribution of p's orig- 
inating from A decay, as a function of the character- 
istics of the original A spectrum. For a A spectrum 
of Boltzmann form with inverse slopes T in the range 
150 MeV < T < 300 MeV, we find that the result- 
ing spectrum of p's within the E917 spectrometer is also 
Boltzmann in shape, with parameters: 



T- 
dN Pde , 



dy 



(0.834 ±0.006)T A 
0.639 *»* 



(0.004 ±0.001) (GeV) 



dy 



The errors reflect the quality of the Boltzmann fits to the 
pdecay spectra, and result in a negligible systematic error 
to the final results. 

The A and p transverse mass spectra were therefore fit 
simultaneously using the functions 



that the A production increases sharply with central- 
ity, from ^ = 1.2lQ gxlO -3 in peripheral collisions to 

^ = 19^5xl0~ 3 in central collisions. The total p pro- 
duction is of comparable magnitude to the A yield, with a 
substantial fraction of the measured p signal attributable 
to decay products from the A. 

We find that the ratio A/p near mid-rapidity increases 
from 0.26lg in peripheral collisions to 3.6^'g in central 
collisions. When combined with the prior observations 
of elevated ratios in central Si+Au collisions [Q and the 
indirect evidence from experiments E864 and E878 Jl5| , 
the present results indicate that A production is larger 
than direct p production in central collisions — and not in 
peripheral collisions — at AGS beam energies. This effect 
is clearly shown in Fig. [jj The mechanism responsible 
for the large A/p ratio is currently unknown. 
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TABLE I: Parameters of an independent fit to the total measured p spectrum (dN Ptotal /dy, Tp total ) and a simultaneous fit (the 
remaining values) to p and A spectra, for the three centrality bins. See text for an explanation of the parameters. Errors are 
statistical followed by systematic. 
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FIG. 1: Invariant spectra as a function of transverse mass for 
A (solid circles) and p (open circles, multiplied by 0.1 for clar- 
ity) in the three centrality classes. The error bars include a 2% 
point-to-point systematic error due to acceptance corrections. 
The solid lines are from the simultaneous fit as described in 
the text and the dashed lines are from independent fits to 
A and p spectra shown for comparison. Solid and dashed fit 
lines overlap for the p spectra. 



3 

>> 
i 

03 

< 



60 



40 



20 



60 



40 



20 



Central 0-12% 




Peripheral 12-77% 




Invariant mass (GeV/c ) 

FIG. 2: The A invariant mass spectra for central and periph- 
eral collisions for the data in Fig. |l[ Solid points are measured 
data, and the histogram is the normalized mixed-event back- 
ground after all corrections have been applied. The A mass 
peak is shaded. 
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FIG. 3: The A/p ratio is shown as a function of the centrality 
given in terms of the percentage of the total interaction cross 
section, a tot . The E917 data are integrated over all p± and 
1.0 < y < 1.4 and are shown as solid circles with la statistical 
error bars for the A/p ratio. The E864/E878 indirect results 
[jn| are at p± ~ and y=1.6 and are shown as most probable 
values with the 98% confidence level of the minimum value 
indicated by the lower error bar. Upper limits are not given. 
Horizontal bars indicate bin widths. 



